Models for Massive Stellar Populations with Rotation 
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O : ABSTRACT 

We present and discuss evolutionary synthesis models for massive stellar populations gener- 
ated with the Starburst99 code in combination with a new set of stellar evolution models ac- 
counting for rotation. The new stellar evolution models were compiled from several data releases 
of the Geneva group and cover heavy-element abundances ranging from twice solar to one fifth 
solar. The evolution models were computed for rotation velocities on the zero-age main-sequence 
> ■ of and 300 km s" 1 and with the latest revision of stellar mass-loss rates. Since the mass cov- 

erage is incomplete, in particular at non-solar chemical composition, our parameter study is still 
preliminary and must be viewed as exploratory. Stellar population properties computed with 
. £^ | Starburst99 and the new evolution models show some marked differences in comparison with 

models obtained using earlier tracks. Since individual stars now tend to be more luminous and 
bluer when on the blue side of the Hertzsprung-Russell diagram, the populations mirror this 
trend. For instance, increases by factors of two or more are found for the light-to-mass ratios at 
ultraviolet to near- infrared wavelengths, as well as for the output of hydrogen ionizing photons. 
If these results are confirmed once the evolution models have matured, recalibrations of certain 
star-formation and initial mass function indicators will be required. 

Subject headings: stars: evolution — galaxies: stellar content — galaxies: individual(SMC, LMC, M 31, 
IC 10) 

1. Introduction still pose serious challenges in evolutionary synthe- 
sis modeling. For instance, evolutionary connec- 
Massive stars are among the most important tions with and between Wolf-Rayet (WR) stars, 
cosmic stellar yardsticks in stellar populations at the si g n i ncance f the nitrogen (WN) and carbon 
low and high redshift. Yet, in many regards they ( WC ) sequence among the WR stars, the ratio 
of blue to red supergiants, stellar surface abun- 

iLaboratoire d'Astrophysique (UMR 5572), Observa- dances , Or stell ar mass loss are all poorly under- 
tone Midi-Pyrenees, 14 Avenue E. Belin, F-31400 Toulouse, , t— — ~ |i „ i J 

-France stood (|Massev 1 12003 ). 
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The Geneva and Padova groups perfected 
the state-of-the-art of stellar evolution model- 
ing in the uppermost Hertzsprung-Russell dia- 
gram (HRD) with the rele ase of complete set of 



Bressan et al. I 19931: Schaerer et al. 1 993allt 



Charbonnel et al. Ill993 



rotation and highlight differences with respect to 
the previous non-rotating grid. We also discuss de- 
tails of the implementation in Starburst99 in this 
section. A parameter study of the most important 
stellar population properties affected by rotation 
is performed in Section 3. In Section 4 we compare 
Bertelli et al.lll994T Fagotto et al.t he model predictions obtained with Starburst99 

to observations. Finally, the conclusions are in 
Section 5. 



stellar evolutio n ary t ra cks (ISchaller et al. _ 1992; 



1994al )bl: iMevnet et al. Ill994l ). These sets have 
become widely used in population synthesis 
modeling (e.g., Fioc fe Rocca-Volmerange 1997t 



2002t Bruzual fc Chariot 20031: Robert et al. 



2003; I Vazquez fc Leitherer 11200 



Leitherer et al. 1999HCerviho et al. 112002 : Schulz et al. 2. Stellar tracks with rotation 



Previous modeling of massive star evolution 
paid particular attention to the size of the convec- 
tive core and stellar mass loss, which were thought 
(together with standard nuclear processing) to be 
the most relevant ingredients in massive star evo- 
lution. This view has changed quite dramatically 
in recent years. It is become clear that stellar ro- 
tation plays a dominant - and in some cases even 
the dominant - role in the evolution of stars wit h 
masses above ~10 M Q (|Maeder fc Mevnet llioooh . 
As a result, the previous generation of stellar evo- 
lution models, which do not take into account ro- 
tation, may be subject to revision. This in turn 
will affect all population synthesis modeling rely- 
ing on these tracks. 

Recently, a new set of tracks has been released 
by the Geneva group (jMevnet fc MaedeT 20051 
and references therein), which addresses stellar 
rotation. This set of tracks should be consid- 
ered "exploratory", i.e., it is still incomplete and 
not meant to fully replace the currently avail- 
able full set of tracks without rotation. These 
new tracks, however, are invaluable in explor- 
ing the fundamental effects caused by stellar ro- 
tation and gauging the impact on stellar popu- 
lation modeling. Consequently we have imple- 
mented this new set of tracks into the latest ver- 
sion 5.1 of the evolutionary synthes is code Star- 
burst99 (jVazauez fc Leitherer 1120051 ). The result- 
ing new version of Starburst99 is not publicly 
available, as the limited number of evolutionary 
tracks does not yet allow the computation of mod- 
els by the astronomical community that could be 
compared to a wide range of observations. 

The present paper is organized as follows. In 
Section 2 we summarize the main astrophysical in- 
gredients in the new stellar evolution models with 



In this section we will describe the most impor- 
tant characteristics of the new set of Geneva tracks 
and highlight differences with respect to previous 
tracks. The discussion will be relatively brief and 
limited to aspects that are directly relevant to the 
implementation of the tracks into Starburst99. 

2.1. Physics of rotation 

The effects of rotation can be classified in three 
categories: 1) the hydrostatic effects, 2) the rota- 
tionally induced instabilities and their impacts on 
the transport of the chemical species and of the 
angular momentum, 3) the effects of rotation on 
mass loss. The hydrostatic effects account for the 
fact that a rotating star is deformed. The stel- 
lar structure equations are modified with respect 
to the non-rotating case . The method used here 
i s the one proposed by iKippenhahn fc Thomas" 



( 197dh . adapted to the case of shel l ular r otation 
as explained in IMevnet fc Maeder I (|l997[) . The 
two main instabilities triggered by rotation are 
the meridional circulation and the shear instabil- 
ities. The effects of these two instabilities on the 
transport of the chemical species and the angu 
lar mo mentum are modeled as proposed by IZahn 
(|l992l ). This theory describes the interactions of 
these two instabilities in a consistent way. It is 
based on one sole a priori hypothesis, namely that 
the star is in a state of shellular rotation (an- 
gular velocity fl constant on isobaric surfaces). 
This hypothesis is quite reasonable from a phys- 
ical point of view, since on an isobaric surface, 
there is no force counteracting the horizontal, i.e. 
isobaric, turbulence. The effects of rotation on 
mass loss have be e n inc luded as prescribed by 
Maeder fc Mevnet ( 2000t ). Note that anisotropic 
stellar winds were included only for the models at 
Z = 0.02. The rotational velocities considered in 
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these models are too modest for this effect to be 
important at other chemical compositions. More 
det ails on the physics of the mode ls can be found 
in (|Mevnet fc Maeder II2003L l2005h . 

2.2. Physical ingredients of the new mod- 
els 

Extensive coverage of the physics of rotation 
in the Geneva tracks is provided in a series 
of papers published by the Geneva group (e.g 



Hirschi. Mevnet fc Maederll2005l : iMevnet fc Maeder 
200.4 and references to prior papers therein) . The 
physical effects of rotation are discussed in detail 
in this series of papers and will not be repeated 
here. The new stellar evolution models used in the 
pre sent population synthesi s mod eling are those 
of (jMevnet fc Maeder lliooll l2005h . The rotating 
models account for the transport of the chemical 
species and of the angular momentum driven by 
meridional circulation and shear instabilities. 

The models were computed with an initial 
equatorial velocity on the zero-age main-sequence 
(ZAMS) of 300 km s _1 . Different processes dur- 
ing the main sequence phase (e.g., the expansion 
of outer layers, transport of angular momentum 
by stellar winds, etc.) decrease this velocity, so 
that the time-averaged velocity on the main se- 
quence for different che mical compositions is 200 - 
250 k m s" 1 (see t ables inlMevnet fc Maeder 1120031 
20051) . Recently, iDufton et al. I (j2006l) have mea- 
sured stellar rotation velocities of Galactic B-type 
stars in NGC 3293 and NGC 4755 whose estimated 
ages are 10 and 15 Myr, respectively. Dufton et 
al. obtained a distribution of velocities that peaks 
at 250 km s^ 1 after accounting for a random dis- 
tribution of the inclination for stars with masses 
between 3 and 1 2 M^. This result is similar to 
that obtain ed bv iHuang fc Gies I (120061). On the 



other hand, Mokicm et al. 



(|2006l ) obtained mean 
velocities of 150 - 180 km s^ 1 for O and B-type 
stars in the SMC with initial masses in the range 
of 15 to 30 M Q . In view of these recent results, 
we see that our averaged velocities might be over- 
estimated for the SMC stars, while they seem to 
be in good agreement for Galactic stars. 

Of course, for a proper comparisons with obser- 
vations, one should compute models for the whole 
range of observed velocities and then convolve the 



results with the distribution of the initial veloci- 
ties. We think that the results of the 300 km s _1 
models will sufficiently well represent the averaged 
evolution so that meaningful comparisons can al- 
ready be performed on the basis of these simple 
models. 

The initial compos i tions were calculated as 



in IMevnet fc Maeder 



( 20011 ). where Y = Y p + 
= 0.23 and AY/AZ = 2.5. 
0.757, 0.744, 0.640 and 



AY/AZ x Z with Y p - 
It follows that X = 
Y = 0.239, 0.248, 0.320 for heavy-element abun- 
dances of Z = 0.004, 0.008, 0.040, respectively. In 
the case of solar chemical composition, the values 
are X = 0.705 and Y = 0.275. The heavy-element 
ratios were assumed to remain constant with Z 
and are the same as those used to calculate the 
opacities for solar composition. The nuclear reac- 
tion rates were taken from the NACRE data base 
( Angulo et al. 19991 ). A more detailed descrip- 
tion of o pacities and nuclear reactions can be also 
found in lVazquez fc Leitherer I (|2005r ). Moderate 
overshooting is include d in all the new models wi th 
and without rotation ( Mevnet fc Maeder 20031 ). 



Fo r the main s e quenc e, the mass-loss rates are 
from lVink et al. I (|2000l 2001), who take into ac- 
count the abrupt change of the ionization and 
opacity of the ste llar wind at certain temper atures, 
as explained by iMaeder fc Mevnet I (l2000l). The 
empirical formula of Ide Jager et al. I ( 1988 ) was 



used f or the domain not covered by Vink et al 
20011 ). For the WR ph ase, the mass-loss rates of 



Nugis fc Lamers I (l2000h were used. We note that 



our adopted main-sequence mass-los s rates were 
derived from theoretical wind models. Puis et al. 



(2006) found good agreement between these rates 
and those observed in O stars in the Galaxy and 
the Magellanic Clouds. The agreement is less sat- 
isfactory, however, for B supergiants whose ob- 
served rates are lower than predicted. These 
rates are empirical, as opposed to the theoreti- 
cal main-sequence rates. The Nugis & Lamers 
rates, which account for the clumping effects in 
the winds, are smaller by a factor of 2 - 3 than 
the rates used in the previous non-rotating "en- 
hanced mass-loss rate" stellar grids. The cor- 
rection factor for enhancing the theoretical mass 
loss due to rotation was taken from Maeder & 
Meynet. The dependence on chemical compo- 



1 On average, v = 4/irvs'mi 



M OC Mq X Z 



0.5 



-WR phases was 


expressed as 


(Kudritzki & Puis 


2000). WR 
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mass-loss rates were assumed to be Z-independent 
in the default set of tracks. However, an addi- 
tional set of tracks with Z-dependent mass-loss 
rates for the WR phase was compu t ed for mod- 



els with Z = 0.040; ICrowther et al. I (|2002a) sug- 



gested the same dependence of the mass-loss rates 
on Z for WR stars as for O stars (oc Z 5 ). The 
mass-loss rates differ from those in the previous 
release of the Geneva tracks. 

2.3. Boundary conditions at the stellar 
surface 

The temperature distribution of a rotating star 
is described by the local radiative flux, which is 
proportional to the effective gravity F oc g e g, 
where g c g is t he sum of the gra vity and the cen- 



trifugal force ( von Zeipel 19241 ). Then the local 



effective temperature on the surface of the ro- 
tating star varies like T e g(v) oc i g c ff(u) 1 / 4 , where 
v is the rotational velocity. Due to this effect, 
the emergent luminosity, colors and spectrum 
will be different for different o rientation angles i 
(jMaeder fc Pevtremann"lll97d ). _n order to avoid 
such complications, we define an average stellar 
T eS by T c 4 ff = L/(aS(Cl)), where a is Stefan's 
constant and S(£l) the total actual stellar surface 
deformed by rotation. This is the average value of 
the effective temperature which will be used for vi- 
sualizing the evolutionary tracks in the theoretical 
HRDs and for evaluating the radiation outputs. 

Since WR winds have a non-negligible optical 
thickness, a correction is applied to new stellar 
models by adopting an effective radius i? e ff at the 
optical thickness r = 2/3. R e g is related to the 
photospheric radius via the relation 



R 



eff 



R 



3k\M\ 

8-KVoo 



(1) 



where n is the opacity, \M\ is the mass-loss rate, 
and is the wind terminal veloc ity. This proce- 
dure is similar to that followed in Schaller et al. I 



([19921 ) in the previous set of tracks. Subsequently, 
the effective temperature at r = 2/3 is calculated 
with the relation L = AttR^ctT^. 

The WR atmospheres are linked to the stellar 
model s using the empirical relation of lSmith et al 
(J2002J) who defined the evolutionary temperatures 
as 0.6T* + 0.4T e ff, where T* is the temperature of 
the WR core. This new definition is a compromise 



between the hotter interior and the cooler surface 
temperature and serves as the interface temper- 
ature used in the expanding model atmospheres. 
Smith et al. empirically determined this relation 
by comparing their model predictions to the ob- 
served temperature distributions of Galactic WC 
and WN stars. 

2.4. The evolution of the rotation velocity 

In the interior of the star, the angular veloc- 
ity varies due to the following effects: contrac- 
tion/expansion of the stellar layers and trans- 
port of angular momentum by convection, merid- 
ional currents and shear instabilities. Without the 
transport processes, the change of f2 would only be 
governed by the local conservation of the angular 
momentum. In that case, there is no coupling be- 
tween the various stellar layers, each one keeping 
its initial angular momentum content. In contrast, 
the transport processes impose some coupling be- 
tween them. During the main-sequence phase and 
in the outer part of the star, the transport mecha- 
nisms relocate angular momentum from the inner 
to the outer regions. This tends to slow down 
the contracting core and to accelerate the expand- 
ing outer envelope. When the stellar winds are 
weak, as for instance in metal-poor regions, the 
angular momentum can accumulate in the outer 
layers so that the surface velocity increases with 
time, approaching the critical limit. When the 
stellar winds are strong, the angular momentum of 
the outer layers is removed and the surface veloc- 
ity decrease, evolving away from the critical limit. 
This illustrates the importance of the mass loss 
for shaping the evolution of the surface velocities 
of massive stars. 

2.5. Consequences for the evolutionary 
tracks 

The differences between the previous and the 
new evolutionary models result from a combina- 
tion of effects. First and foremost, rotation in- 
troduces the majority of the changes. In addi- 
tion, the opacities and overshooting are different. 
The former produce a shift in temperature towards 
lower effective temperatures. The overshooting is 
responsible for lower effective temperatures dur- 
ing the main-sequence, and combined with a core 
extension due to rotation it also affects the blue 
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loops during the He-burning phase. Other major 
effects are: 

• Tracks with rotation become more luminous 
than those with no rotation. As a conse- 
quence, the derived mass of a star on the 
HRD will be lower than compared with non- 
rotating tracks. 

• The main-sequence width is enlarged due to 
rotation and overshooting for stars with ini- 
tial masses below ~ 30 M©. This effect re- 
sults from a larger He-core produced during 
the H-burning phase. 

• For masses higher than 50 M Q , the opposite 
effect occurs. For instance, the reddest point 
reached by the 60 M non-rotating stellar 
model at the end of the main-sequence phase 
is log T c ff ~ 4.15, whereas for the rotating 
model, it is ~ 4.6. This is a consequence of 
the lower diffusion time scale ratio for chem- 
ical elements at higher initial mass. Mas- 
sive stars will be more mixed than stars with 
lower masses at the same evolutionary stage. 
A higher degree of homogeneity is reached, 
and this turns the star bluer. Combined 
with these effects, rotation and mass loss 
change the stellar surface composition and 
favor the formation of WR stars. 

• The mass-loss rates used in the present mod- 
els are smaller than t he high rates used in 



Mevnet et al. I (|1994[ 1. Therefore, the fi- 



nal masses at the end of the He-burning 



phase are larger than those in lSchaller et al. 
(Il992h . 



• The main-sequence life times at Z = 0.020 
and 0.04 have increased by 15 - 25 and 16 - 
27% in comparison with the non-rotating 
models, respectively. This is mainly an ef- 
fect of rotational mixing which continuously 
supplies the core with fresh nuclear fuel, i.e., 
hydrogen. The increase of the luminosity 
due to rotation cannot compensate for this 
effect, and the main-sequence life times are 
increased by rotation, at least in this metal- 
licity range. 

2.6. Predicted properties of WR stars 

The WR phase is significant for massive stars, 
and stellar rotation can dramatically affect this 



phase. Our knowledge of this phase is in- 
sufficient to permit constraints on the proper- 
ties of WR stars purely from first principles 
(|Mevnet fc Maeder |[200a 120051) . Therefore some 
ad hoc assumptions must be made. For the present 
stellar models, the star enters into the WR phase 
when 1) logTtff > 4.0 and 2) when the mass frac- 
tion of hydrogen at the surface is X s < 0.4. There 
is no significantly different outcome if the value for 
X s varies bet ween 0.3 or 0.4. In the first release of 
Starburst99 ( Leitherer et al. Il"l999l ). the definition 



of a WR star was based on the work of Maeder 



(|199ll ) , who used a minimum temperature limit for 
a WR star of logT cff > 4.4. This limit is 0.4 - 0.5 
dex higher than that assumed in the new tracks. 
Recent observations of WC/WN ratios in the Mag- 
ellanic Clouds resulted in b etter agreemen t with 
a lower temperature limit ( Massey 1 120031 ) . The 
WR definition via Teff is arguably rather adhoc 
and needs to be addressed. Therefore we have 
studied synthesis models with both temperature 
limits and with the old and new sets of tracks to 
gauge the importance of this parameter for the 
WR properties. 

We are also interested in the stellar mass at 
the entrance and during the WR phase. As a 
reminder, the most massive stars enter the WR 
phase while still in the H-burning phase. Two dif- 
ferent mass limits for the W R phase have been 



define d by t he Geneva group iMevnet fc Maeder 



(|2003l l2005l) . The first one is the mass limit 



to enter to the WR phase on the main-sequence 
(Mowr); the second mass limit is the mass to en- 
ter the WR phase at any time during the evolu- 
tion of the star (Mwr). The first limit implies 
an evolutionary scenario avoiding the luminous 
blue variable (LBV) phase so that the phases cov- 
ered are O-eWNL-eWNE-WC/WOH The second 
mass limit implies an evolutionary scenario first 
going through the O-LBV or red supergiant (RSG) 
phases and then following the sequence cWNL- 
eWNE-WC/WO. The mass limits for different sets 
of tracks are shown in Table [TJ No rigorous com- 
parison with observations has been made but the 
sparse data suggest a h igher Mowr clo ser to 75 
Mq for Galactic stars (jCrowther 



2007). Fig. CD 



2 The transition from the WN to the WC phase occurs when 
the mass fraction of carbon is higher than 10% of the mass 
fraction of nitrogen, "e" denotes a spectral su bclass defined 
by evolutionary criteria l|Foellmi et al. II2003T ). 
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shows evolutionary tracks for 60 M Q with an in- 
dication of the entry point into the WR phase for 
four different chemical compositions for the new 



models with u rot 
models with v rot 



= 300 km s" 



Thc tracks for 



= km s 1 and the previous 



tracks with high mass-loss rates are also included 
in Fig. [lj 

2.7. Evolutionary phases and implementa- 
tion into Starburst99 

The principal new set of tracks has an initial 
i> rot of 300 km s _1 with four different values of the 
chemical composition, as described in Section 2.2. 
A second set of tracks has v ro t — km s _1 . The 
non-rotating models were used to calibrate and 
compare the new models with rotation to the pre- 
vious sets of tracks. These tracks were generated 
for Z = 0.02 and 0.04. An additional set of tracks 
with Z = 0.04 and metallicity dependent WR 
mass-loss rates is used in this paper as well. 

The tracks were computed up to the end of the 
He-burning phase in all cases. Table [2] summarizes 
the covered mass range for the different metal- 
licities. The tracks with solar chemical composi- 
tion have the most complete coverage of parame- 
ter space. Tracks were released in the same way as 
the previous set of tracks by considering equivalent 
evolutionary time-step points along each track. 
The number of points, however, increases by a fac- 
tor of 7. Point 1 defines the ZAMS, point iOO is 
the end of the H-burning phase, point 20f is the 
beginning of the He-burning phase, and point 350 
is the end of the He-burning phase. 

The new set of tracks predict the same stellar 
parameters as in the original set of the Geneva 
tracks. New parameters like the ratio of the po- 
lar to the equatorial radii, equatorial velocity and 
angular velocity (ft/ftbreakup) are also included in 
the new release. Throughout the paper, we will 
use the following nomenclature: the "old Geneva 
tracks" are the previous tracks with high mass-loss 
rates used as default by Starburst99. The "new 
Geneva tracks with or without rotation" are the 
new tracks with v ro t — 300 km s _1 and km s _1 , 
respectively. 

The method of implementation into Star- 

I u 

burst9 9 was similar to the one used bv lVazquez fe Leitherer 

(2005). However, the number of grid points of the 

new tracks is different. In order to be consistent 



with the existing, old set of tracks, we carefully 
selected 77 grid points to define key evolutionary 
phases and cover periods of rapid evolutionary 
changes. To achieve consistency, we found it nec- 
essary to resample the old Geneva tracks from the 
original 51 to 84 grid points. The old tracks ex- 
tend up to the C-burning phase. We re-built the 
new tracks with 84 points as well, even though 
their evolution does not go beyond the He-burning 
phase. This leaves 7 slots available for a future 
implementation of the C-burning phase. 

3. Rotation versus non-rotation 

The purpose of this section is to compare inte- 
grated parameters from Starburst99 obtained with 
the old and new tracks. Constant star forma- 
tion and instantaneous burst models were gener- 
ated with an initial mass function (IMF) index of 
a = 2.35 over the mass range of 1.0 - 100 M 
until 100 Myr. We used the fully line-blanketed 
atmospheres of ISmith et al~l (l2002h for hot stars 



and A TLAS9 models adjusted by iLeieune et al 



( 1997L 1998) for stars of spectral type B and later 



3.1. Isochrones 
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As a first comparison, we show in Fig 
isochrones for two different abundances (Z = 
0.008 and 0.020) produced with the new tracks 
with rotation of v YO t — 300 km s _1 (left) and those 
produced with the old Geneva tracks (right). In 
each case the oldest isochrone corresponds to the 
age at which the lowest mass included in the new 
models with v Iot — 300 km s _1 is present (see 
Table EJ. 

Isochrones produced with the new tracks are 
more luminous than those obtained with the old 
Geneva tracks. Stars on the RSG branch with 
Z = 0.008 are bluer than those obtained with the 
old Geneva tracks. In the case of Z = 0.02, the 
new isochrones do not have the extended blue loop 
except for the model with 9 M Q . As this is the low- 
est mass model included in the new set of tracks, 
this isolated blue loop for just one mass will cre- 
ate a bump in the integrated light pr o pertie s, as 

(|2003h dis- 



we will see below. Mevnet fc Maeder 



cuss how the reduction of the blue loops in the new 
acks with rotation results from the combined ef- 



fect of overshooting and the core extension by ro- 
tation. 
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3.2. Spectral energy distributions 

We have compared the stellar spectral energy 
distributions (SEDs) for young ages. Throughout 
this paper we follow the convention to plot the 
models with the different sets of tracks from the 
Geneva group with enhanced mass loss, the new 
models with v rot — 0, and the new models with 
D rot = 300 km s _1 with solid, short-dashed and 
dotted lines, respectively. Young population mod- 
els with Z — 0.020 in the range of 1-15 Myr were 
divided and are plotted in two different figures. 
Fig'-Elshows models up to 6 Myr, and Fig.|4]shows 
the models with older ages. SEDs produced by 
the models with rotation show an energy increase 
around the neutral helium edge with respect to the 
models with old Geneva tracks, which can be up 
to ~ 5 dex. This is the result of a higher effective 
temperature for a given age in the new models. 
The flat sections at the shortest wavelengths of 
some spectra in Fig. [3] and Fig. 2] are numerical 
artifacts. 

The difference between the old and new models 
in the range of ~ 2.4 < log A < 2.95 at ages 3 - 
7 Myr may not appear very striking with the scal- 
ing in Figures [3] and [U Nevertheless, it amounts 
to about a factor of two for the number of ionizing 
photons, as we will discuss below. The onset and 
disappearance of the enhanced ionizing luminos- 
ity coincides with the presence of WR stars in the 
new models with rotation. 

At even longer wavelengths above ~ 10000 A 
(Figure [5]), the difference in luminosity is < 0.5 
dex for the same young ages covered in Figs. [3] 
and SI When the age increases and RSGs begin 
to dominate the red and near-infrared part of the 
SED, the models with rotation become even more 
luminous as they predict a larger number of RSGs. 

We restrict our discussion of the SEDs to the 
case of Z = 0.02 because only this set of evolu- 
tionary models is sufficiently complete to allow the 
computation of a meaningful SED. 

3.3. Massive-star inventory 

3.3.1. The WR/O ratio 

The first test for studying the effects of rotation, 
mass loss, and temperature change is the ratio of 
the number of WR to O stars for an instantaneous 
burst. Fig. [5] shows this ratio for six different in- 



stantaneous burst models: old Geneva tracks and 
log T c ff > 4.4 (red solid line); old Geneva tracks 
and logT c ff > 4.0 (red long-dashed); new tracks of 
Vmt = 300 km s" 1 and logT ff > 4.4 (blue dotted); 
new tracks of v m t = 300 km s _1 and logT e ff > 4.0 
(blue dot-dashed); new tracks of v m t = km s _1 
and logTcff > 4.0 (black short-dashed); new tracks 
of v TO t = 300 km s" 1 , logT e ff > 4.0 and different 
mass mass loss rate for WR phases at Z = 0.040 
(dot long-dashed). As a word of caution, the new 
models with lower chemical composition than so- 
lar are incomplete at the lowest masses. While the 
WR evolution is well covered using the new tracks, 
the evolution of the least massive O stars is not. 

The WR-rich phase after an instantaneous star- 
burst covers a larger interval of ages when rotation 
is accounted for. This comes from the fact that 
rotational mixing aids an entry at an earlier evo- 
lutionary stage. In order to appreciate this point 
better, recall that mass loss by stellar winds is the 
key process allowing nuclearly processed layers to 
appear at the surface and thus make the star en- 
ter into the WR phase in non-rotating models. In 
the case of rotating models, another effect con- 
tributes in building the chemical characteristics of 
the stellar surface of WR stars: the mixing in- 
duced by rotation brings elements produced in the 
stellar core to the surface. Thus, for the present 
mass-loss rates, rotation favors an entrance into 
the WR regime at an evolutionary stage when the 
central hydrogen mass fraction is higher. As a con- 
sequence, the WR lifetimes will be longer. Also, 
stars with lower initial mass can be progenitors 
of WR stars. An immediate consequence are the 
peaks in the WR/O ratio between 6 and 9 Myr 
for Z = 0.020 and 0.040. These artifacts result 
when WR stars exist but most O stars have al- 
ready evolved into B stars. 

Furthermore, the WR phase lasts significantly 
longer in the new models so that the overall effect 
is a net increase of the WR/O ratio in an equilib- 
rium population. This becomes obvious in Fig. [71 
where we have plotted the same ratio for constant 
star formation. The higher number of WR stars 
in the new models is significant. 

Generally, there are more WR stars in the new 
stellar models with rotation and with the lower 
temperature limit for the WR phase. It is impor- 
tant to realize that most of this difference is caused 
by WR stars evolving from the lowest masses. 
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These are the most populous species, and varying 
their numbers has the most dramatic effect. 

3.3.2. The WC/WN ratio 

The same cautionary remarks on incomplete- 
ness as in the previous subsection apply here. 
Fig. [5] shows the WC/WN ratio for instantaneous 
star formation. As in the case of the WR/O ratio, 
we plot models for constant star formation as well 
(Fig. [9]). The general trend in these figures is a 
decrease of WC/WN, in particular at low chemi- 
cal composition. This effect is most dramatic for 
models with Z = 0.008 and 0.004. As discussed by 
Massev I ([20031 ) , the previous generation of Geneva 
tracks has been plagued by predicting a too high 
WC/WN ratio for metal-poor environments. The 
new models appear to be an improvement in this 
respect (cf. Fig. EH). 

The dot long-dashed line plotted for Z = 0.040 
in Figs. [8] and [9] is a model for the WR phase with 
the new tracks of u ro t = 300 km s _1 using the 
lower temperature limit (log T c ff > 4.0) and a dif- 
ferent prescription for the mass-loss rate during 
the WR phase. These tracks were specifically in- 
troduced to boost the predicted number of WC 
stars for higher chemical composition. However, 
contrary to expection, Figs. [8] and [9] suggest that 
these models do not change significantly the re- 
sult: the WC/WN ratio is slightly lower (dot-long 
dashed line) than in the corresponding new models 
without Z-dependent WR mass-loss rates (dotted 
line). 

3.3.3. The R/WR ratio 

Another often employed ratio is that of RSGs 
over WR stars (R/WR). Before exploring this ra- 
tio, it needs to be defined first. O ur definition is 
based on the work of lMassev I (|2002h who considers 
as RSGs all stars with (V— R)o > 0.6, correspond- 
ing to stars of logT cff < 3.66 (4600 K, or late G 
type) and at the same time having Mboi < —7.5. 
This definition ensures that only very luminous 
red stars are counted and less luminous AGB stars 
are excluded. We have implemented this definition 
in Starburst99 to calculate the predicted number 
of RSGs. 

As in the previous sections we are taking into 
account different temperature limits for the WR 
phase. We do not discuss models with instanta- 



neous star formation as R/WR is indeterminate in 
this case. The WR and RSG phase are not over- 
lapping in the evolutionary models, and the ratio 
of the two quantities would be either zero or infin- 
ity. This degeneracy results from our extreme as- 
sumption of an idealized single stellar population. 
In the presence of a small age spread at star birth, 
there is a small likelihood of both RSGs and WR 
stars co-existing. Observations of clusters suggest 
few if any such examples. Fig. [TUl shows the R/WR 
ratio for constant star formation. To understand 
this figure, we have to take into account the life 
times of massive stars, which enter when stellar 
populations in equilibrium are considered. 

The panel for Z — 0.020 suggests values of order 
unity for R/WR, with only a mild dependence on 
the adopted tracks or temperature definition. In 
the case of Z = 0.040, the spread between the dif- 
ferent tracks becomes much larger. The difference 
can be up to ~ 1 dex, independently of the temper- 
ature limit for WR formation. The new tracks for 
this chemical composition predict a higher number 
WR stars (therefore a lower R/WR ratio) because 
the WR phase sets in very early close to the main 
sequence. 

The panels for Z = 0.004 and 0.008 should be 
interpreted with care because of the incomplete- 
ness of the tracks. The differences between the 
individual curves are caused by variations of the 
WR numbers only. Massive stars with M > 40 
M© will never reach the RSG phase in this part 
of the HRD. We will return to this point when we 
analyze the HRD and compare it with observa- 
tions. On the other hand, those RSGs that would 
contribute to R/WR are mostly stars with initial 
masses < 25 M Q but they are not included in the 
new models because of incomplete tracks (cf. Ta- 
ble E). 

3.3.4. The B/R ratio 

Next, we discuss the blue-to- red supergiant ra- 
tio (B/R). As in the previous section, we need to 
define th is ratio first. W e will again use the defi- 
nition of iMassev"! ( 20021 ). who defines blue super- 
giants as having (B — V)o < 0.14, corresponding 
to stars of logT cff > 3.9 (8000 K, or type A9 I), 
and Mboi < —7.5. As in the case of RSGs, we im- 
plemented this definition into Starburst99 to cal- 
culate the number of blue supergiants. 
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Fig. ITTI shows the B/R ratio for constant star 
formation. The case of instantaneous star forma- 
tion has been omitted for a similar reason as be- 
fore: very luminous blue supergiants have no evo- 
lutionary counterpart in the uppermost HRD. 

The trends in the individual panels are qualita- 
tively similar to those for R/WR. The two panels 
for Z = 0.004 and 0.008 are again to be taken 
with care. The apparent increase of B/R with de- 
creasing Z is an artifact of incomplete tracks. In 
contrast, the panels for Z = 0.020 and 0.040 re- 
flect the action of the new evolution models and 
are suitable for comparison with observations. 

3.4. Ionizing continuum 

Figure [12] shows the number of ionizing photons 
below 912 A, 504 A, and 228 A for an instanta- 
neous starburst. The photon output of the three 
sets of tracks is similar during epochs when WR 
stars are not present. In contrast, significant dif- 
ferences are observed in the presence of WR stars 
between 3 and 8 Myr. The models with rotation 
generate more ionizing photons than the original 
tracks. The effect is most pronounced for the ion- 
ized He continuum, which was opaque with the old 
tracks but becomes more luminous by many orders 
of magnitude in the new models. Note, however, 
that the photon output in absolute terms is still 
very small for most of the time and may not be 
observable. The effect does become observable in 
the neutral He continuum, as well as in the neutral 
hydrogen continuum. The latter displays a factor 
of 2 - 3 increase over the old models. We empha- 
size the immediate consequences on, e.g., the cal- 
ibration between the ionizing luminosity and the 
star-formation rate and its dependence on Z, as 
the WR/O ratio is strongly Z dependent. 

As a reminder, the treatment of T e g at the in- 
terfa ce with the extended a tmospheres is the same 
as in ISchaller et al. ] (|l992h . but i? ft at the loca- 
tion of the optical depth r(R c ff) — 2/3 changes 
according to the new rotation physics. Further- 
more, we have modified the temperature used 
in the expanding a tmospheres as suggested by 
Smith et al. 1 |2002). Using a different recipe for 
the interface temperature would change the prop- 
erties of the WR population. Although the recipe 
for the temperature is the same as in our previ- 
ous Starburst99 model set, the WR parameters of 
the evolution models are quite different now. This 



explains why the ionizing photon output is very 
different, despite using the same WR model at- 
mospheres. 

3.5. Panchromatic luminosities 

One of the most important consequences of ro- 
tation is that the new stellar models are more lu- 
minous at any wavelength. We compared the pre- 
dicted luminosities from an instantaneous burst 
for different spectral bands. As before, the old 
Geneva tracks, new tracks of v Iot = 300 km 
s _1 , and t> rot = km s -1 are represented by 
solid, dotted, and short-dashed lines, respectively. 
Figure [TBI shows the bolometric (M bo i), 1500 A 
(Li 500 ), B-band (Mb), V-b&nd (My), J-band 
(Mj), and JT-band (Mk) luminosities for solar 
composition. The photometric system, the L1500 
definition and the zero poin t for Mboi are as in 
Vazquez fc Leitherer I (|2005l) . 

We see from Fig. [13] that the luminosities are 
more affected by the new models with i> rot = 300 
km s _1 than by those with w rot = km s _1 . 
The contribution to the bolometric and short- 
wavelength bands is mainly due to hot main- 
sequence and WR stars for young populations with 
ages < 15 Myr. In addition, all tracks have higher 
bolometric luminosity, thereby raising Lisoo- At 
longer wavelengths, the increased emission is due 
to the excursion of stars with masses < 30 M Q 
to the RSG phase for ages > 20 Myr. The RSG 
phase is significantly enhanced in the new models 
with rotation. 

3.6. Colors 

Fig. [JJ] displays the color evolution for an in- 
stantaneous burst at three different chemical com- 
positions. The usual caveat as to incompleteness 
at low Z applies. We decided to include results for 
metal-poor populations as they can provide im- 
portant astrophysical insight when interpreted in 
context. The readers are urged to keep in mind 
the compilation of tracks in Table [5] 

The contribution from RSGs to the red bands 
is most obvious in (V — R) and (V — K) from 
models with w rot = 300 km s _1 after 15 Myr. The 
effects on colors are delayed with respect of the 
old models. The new models with i; rot = km 
s _1 follow the prediction of the old Geneva tracks 
except for younger ages < 10 Myr, when the new 
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tracks evolve faster than the old ones. 

In general, for all plots shown in Fig. [HI the 
new models with w rot = 300 km s _1 predict values 
of (U — B) bluer by up to ~ 0.2 magnitudes than 
the models using the old Geneva tracks. There is 
a higher contribution in the U band for popula- 
tions < 7 Myr. The (B — V) color is less affected, 
but still shows a noticeable effect. The difference 
between old and new models can be up to ~ 0.1 
magnitudes during some epochs. 

The major differences at solar chemical com- 
position occur for (V — R) and (V — K). After 
10 - 15 Myr these colors are redder by up to 1.5 
magnitudes in the new models with rotation. As 
discussed before, this effect is due to the redward 
evolution to the RSG phase of stars with masses 
lower than 30 Mq. In the new models, the stars 
spend a longer fraction of their life times in the 
red part of the HRD, rather than evolving back 
to higher temperatures on blue loops. The peak 
seen around log(Age[Myr]) ~ 1.6 in Fig. [Mb is 
the effect of the blue loop for the model of 9 M Q . 
This loop is introduced by the interpolation in the 
mass range 7 < M/M Q < 12 (see the discussion 
in Section 2). 

3.7. Mass loss and light-to-mass ratio 

In the absence of rotation, the prescription used 
to predict the mass-loss rate leads to lower rates on 
the main-sequence than in the old Geneva tracks. 
In combination with rotation, the mass-loss rates 
are enhanced by the time the stars enter the WR 
phase at the end of the main-sequence. There- 
fore, the total mass loss and all related parame- 
ters, such as the mass-to-light ratio, are different 
in comparison with the previous models. We have 
plotted in Fig. Qjj] the mass-loss rate of instanta- 
neous burst models using the new tracks with ro- 
tation of v Ia t = 300, new tracks of v rot — km 
s _1 and models using the old Geneva tracks la- 
beled as dotted, short-dashed, and solid lines, re- 
spectively. We have also included the set of tracks 
with Z = 0.040 with a mass-loss rate oc Z 5 for 
the WR phase. This set is plotted as the dot- 
dashed line in the bottom right panel labeled as 
Z = 0.040. 

Fig. [15] suggests somewhat higher mass-loss 
rates close to the ZAMS for models with rotation 
than for models with the old Geneva tracks. This 



effect is more pronounced at lower chemical com- 
position. Subsequently, the new mass-loss rates 
are lower at the end of the main-sequence but will 
again surpass the rates of the original models in 
the WR phase. It is important to recall that the 
comparison in Fig. [15] is for the mass-loss rates in- 
tegrated over the entire population at a given evo- 
lutionary phase, and not for the mass-loss rates 
at a given position in the HRD. The rates used 
for the new models with rotation during the WR 
phase would be lower than those used in the old 
Geneva tracks for a given position in the HRD. 
In the case of the new rotating models, since the 
stars have higher mass when they enter the WR 
phase, they have also higher luminosities and thus 
higher mass-loss rates. 

The figure for Z = 0.040 includes an additional 
model (dot-dashed line). This model is based on 
tracks using a mass-loss rate oc Z 5 for the WR 
phase. This scaling produces somewhat higher 
rates in the WR phase at super-solar chemical 
composition, an effect that can be clearly seen in 
the figure. While this has important consequences 
for the evolution to and within the WR phase, the 
impact on the overall mass return is small. 

The total mass loss is shown in Fig. [16] This is 
simply a time- integrated version of Fig. 1151 and its 
behavior mirrors that of Fig. Qjj] The line types 
are the same as in the previous figure. Despite the 
differences between the mass-loss rates early and 
late in the evolution, the total mass lost over the 
entire stellar life time is rather similar in all sets 
of evolution models. 

Fig. fT7l shows the light-to-mass ratio for Lbolj 
Ly, Lk- Panels a) and b) are for Z = 0.020 
and Z = 0.040, respectively. Line types are as 
in Fig. [16] The model prediction for the bolomet- 
ric light-to-mass ratio is higher by ~ 0.2 dex in the 
models with v mt = 300 km s~ x for both chemical 
compositions once post-main-sequence stars dom- 
inate the light. The fact that there is no differ- 
ence between the models with Z = 0.020 after 
log(Age[Myr]) = 1.7 is an artifact. This age cor- 
responds to the life time of a 9 M star, which 
equals the lowest-mass track available for the new 
evolution models. After this age, the tracks for all 
population models are the same. At Z = 0.040, 
the same effect sets in after about 10 Myr. The 
increase of the L/M ratio results from the higher 
luminosities of the new models with rotation. As 
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we showed before, the total mass loss, and there- 
fore the current stellar masses are very similar in 
the model sets. The increase of the L/M ratio 
during the post-main-sequence phase in the new 
models is even more pronounced in the V, and in 
particular in the K band. The passbands favor 
the contribution by luminous post-main-sequence 
stars, whereas the bolometric luminosity tends to 
have a strong component from hot main-sequence 
stars whose luminosities are not too different in 
the three model sets. 

L/M for the bolometric luminosity with tracks 
of frot = km s _1 and with the old tracks are 
essentially identical. Some differences can be rec- 
ognized for the light-to-mass ratio in the V and 
K bands. There is a higher contribution around 
logAge[Myr] = 0.7 (Fig. PTTb). In general, the 
non-rotating tracks evolve faster and are less lu- 
minous than models of v rot = 300 km s _1 . There- 
fore, the feature is caused by the earlier excursion 
to the red part of the HRD in these models. 

Fig. [TTb includes the prediction of the models 
with tracks using a mass-loss rate oc Z°- 5 (dot- 
long dashed line). These tracks are sampled for 
masses in the range of 40 < M < 120 M©. As a 
result, their contribution disappears after 10 Myr. 
The models produce significantly higher luminosi- 
ties when WR stars are present and contribute to 
the light. This is the case in particular for the V 
band. 

4. Comparison with observations 

Having compared the new models to the pre- 
vious generation of non-rotating models, we are 
turning to a first comparison with observations. 
The incompleteness of the new grid limits this 
comparison to relatively few observables. Never- 
theless, these early indications can provide guid- 
ance for future modeling efforts. 

4.1. Comparison with synthetic Hess dia- 
grams 

Hess diagrams are useful tools for visualizing 
actual number densities in HRDs, as opposed to 
simply plotting evolutionary tracks. This is par- 
ticularly important for the comparison of stellar 
life times. Ordinarily, stellar number densities 
are represented in Hess diagrams. In the follow- 
ing we have generalized this method by plotting 



integrated population properties such as rotation 
velocities or abundances weighted by the stellar 
number densities. 

The Hess diagrams were calculated with a 
grid defined by 56 points over the range 4.9 > 
log TV/ > 3.5 and 2.5 < log L/L Q < 6.5. Star 
formation was assumed to be continuous over 
100 Myr, which implies star formation and stel- 
lar death are in equilibrium for masses as low as 
~5 M Q . A Salpeter IMF in the range of 1 - 100 
Mq was used in combination with the new tracks 
of v TO t = 300 km s _1 . Variables visualized in these 
diagrams are: star number (iV stars ), current mass 
(M c ), mass- loss rate (log dm/dt), mass fraction of 
H, He, C, N, O, and the rotation velocity. The 
diagrams were made with an IDL routine that 
uses the subrouti ne sinqlecontour ' from the pack- 
age CHORIZOS (|Maiz-Apellamzll2004h . 

Stellar statistics from star-forming galaxies in 
the Local Group such as the Magellanic Clouds, 
are well suited for comparison with ou r mod - 
els. We have used the results of Venn (1999), 



Venn fc Przvbilla I (l2q03fl.lCrowther et al I(l2002bft. 
fouret et al. I (l2003h.lHillier et al. l(l2003flJLennon et al. 



2003ft , ITrundle et al. I (|2004ft . and lTrundle k. Lennon 
(|2005l ) who analyzed samples of massive stars in 
the Small Magellanic Cloud (SMC). The data, 
including the errors, are shown in Table [3l We 
compared their nitrogen abundance, current mass, 
mass-loss rate, and rotation velocity to those pa- 
rameters in the new Geneva tracks of v Iot — 300 
km a~ 1 via the Hess diagrams. 

Fig. HU shows the Hcss-diagram at Z = 
0.004 for the nitrogen abundance in units of 
log[(iV/iJ)/(Ay_ff)mi]. For a fair comparison with 
the data we have re-scaled the nitrogen abundance 
from the tracks to the average a bundance of ni 
trogen of the SMC as done by ITrundle et al. 
(|2004ft . The resulting ZAMS abundance is 
12 + log(iV/iJ)mi = 6.55. The grey-scale visu- 
alizes the average nitrogen abundance in each 
temperature-luminosity bin predicted by the mod- 
els under the assumption of constant star forma- 
tion. The average value can be the result of the 
superposition of different evolutionary phases. For 
instance, an evolutionary track may cross the same 
bin on its right- and leftward movement, with the 
star having different chemical composition. In or- 
der to facilitate the comparison with the data, we 
chose matching colors for the symbols represent- 
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ing the data. Each symbol is surrounded by a 
square. In those cases where the models and data 
agree very well, the symbol is invisible and only 
the square appears. 

In the case of Z = 0.004 (as well as Z = 0.008), 
the new tracks cover the uppermost mass range 
only. We substituted the old Geneva tracks for 
the missing rotating tracks in order to be able 
to generate the Hess diagram. (This was done 
in Fig. [TB] and all other figures with missing mod- 
els.) The cut-off luminosity below which the tracks 
were substituted is indicated by the horizontal 
line. The data added to Fig. [18] are from Ta- 
ble [3] The average nitrogen abundance for all 
stars in this table except the A and F supergiants 
is 12 + log N/H ~ 7.64. Very few stars are lumi- 
nous enough for a comparison with the new mod- 
els. Taken at face value, those few stars suggest 
the nitrogen surface abundance is somewhat un- 
derestimated by the new stellar models. 

In Fig.[19]we show the Hess diagram for the cur- 
rent stellar mass. The "observed" masses were de- 



those predicted from the models. The average ob- 
served value for massive stars (ignoring A and F 
supergiants) is usini ~ 65 ~ km s . 



termin ed spectroscopically by lTrundle fe Lennon 
(|2005l ). Note that these masses were derived 



from an atmospheric analysis and not dynamically 
from binary stars. Therefore the term "observed" 
should be taken with care. The data agree very 
well with the model predictions, with the caveat 
that the comparison for most of the stars refers to 
the old non-rotating tracks. On the basis of the 
limited sample we conclude there is no discrepancy 
between evolutionary and spectroscopic masses. 

The Hess diagram for the mass-loss rates is in 
Fig. [20] The mass-loss rates are of course closely 
related to the stellar current mass (Fig. Q1J]) as 
these rates determine the relation between the 
masses at any epoch during stellar evolution. The 
new models match the observation very well in 
most stars over the range of — 7 < log dM/dt < 
—6, with a few stars below this value and some 
above. Most notably, the predicted rates for B 
supergiants are higher than observed. 

Next we turn to the rotation velocity. Fig. I2T1 
shows the Hess diagram for the rotation velocity 
v rot for the stellar sample in Table [3] The data 
were corrected by the statistical factor for the dis- 
tribution of the inclination angle 4/7r to obtain 
?; ro t . Therefore the values of the grey-scale are the 
rotation velocity at the stellar equator. Fig. [21] 
suggests lower observed rotation velocities than 



Trundle et al. I (|2004h compared their derived 



rotation veloci ties with an early se t of tra cks with 



rotation from iMaeder fe Mevnet I (|2001 1 . In or- 
der to enforce agreement, Trundel et al. had to 
postulate an average inclination angle of i = 20°, 
which means most stars are observed pole-on. Al- 
ternatively the observations could be systemati- 
cally low, or the models too high. Apparently, the 
discrepancy remains with the new set of tracks, 
as suggested by Fig. [21] Trundle et al. discuss 
a possible overestimate of the rotation velocity in 
their sample because the line widths could be af- 
fected by macroturbulence, which would add to 
the broadening. However, it seems doubtful that 
macroturbulence can be sufficiently large to ac- 
count for the discrepancy. A more likely explana- 
tion could be a selection effect: the star s in Trun- 
dle et al. and lTrundle fc Lennon I ([20051 ) were cho- 
sen to have low rotation velocity (Lennon, private 
communication) . 

The comparison of the rotation velocities with 
SMC stars suffers from the incompleteness of the 
evolutionary models. Since the models with so- 
lar chemical composition are complete down to 
9 M Q (Table 2), we performed a corresponding 
comparison for t h e sam ple of Galactic O stars of 
Conti fc Ebbetsl (|l977l ). We translated the spec- 



tral ty_p^s_to_tem£exatm^and luminosity by using 
the Schmidt-Kaler I |l982) calibration. Spectral 
type bins occupied by more than one star were 
averaged in vsini. The 119 stars from the sample 
produced 33 averaged values in the Hess diagram 
shown in Fig. [55] The values are also tabulated in 
Table [4] Once corrected by 4/7r, the observed val- 
ues remain somewhat smaller than the model pre- 
dictions but the difference is much smaller than 
that for the SMC stars in Fig. [21] Stars with 
masses lower than 30 M show a larger discrep- 
ancy. 

We also compared the rotation velocities of 
the solar abundance m odels to those derived by 
iHuang fc Gies (2006). Their sample consists 
mostly of B stars in 19 young open clusters, which 
is a useful extension of the previous O-star sam- 
ple to lower luminosities. As before, we averaged 
the rotation velocities of identical spectral types. 
These averaged values are shown in Table [5] and 
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plotted in Fig. [23j Once again, we have corrected 
the observations by the statistical factor of 4/ir to 
compare with the models. As for the O stars, the 
models somewhat overestimate the rotation veloc- 
ities in this collection of Galactic B-type stars. 

To summarize, we find good agreement between 
the models with rotation and the data for the cur- 
rent masses and the mass-loss rates. The agree- 
ment for the nitrogen abundances and rotation 
velocities is poor but the observational material 
is sparse and subject to selection effects. 

4.2. Red supergiants in the Magellanic 
Clouds 

One of the main drivers for building a new set 
of tracks was the long-standing problem of the 
disagreement between t he predicted and observed 



RSG num ber statistics (Langer fe Maeder 1995; 



Maeder fe Mevnet 11200 II ). In this section we com- 



pare the predictions of the new set of tracks with 
rotation to samples of RSGs in the Magellanic 
Clouds. 

(|2003l) 



Massev fe Olsen 



conducted a photo- 
metric survey of the Magellanic Clouds using 
BVR CCD photometry and found 105 and 158 
confirmed RSGs in the SMC and LMC, respec- 
tively. These authors performed a comparison 
of their census with different evolutionary models 
and found that none of the tested models matched 
the RSG sequence. There are two issues involved 
in this apparent failure: first, the observed and 
predicted effective temperatures do not match. 
This mismatch may be mitigated or even com- 
pletely resolved by adopting a revised, hotter T e g 
scale. A second, more troublesome issue are the 
lifetimes. The time spent in the very red part of 
the HRD is too short in the models. As a result, 
the stellar number densities in the region where 
RSGs are observed do not match. However, as 
pointed out before, most RSGs are produced by 
stars with masses < 30 M Q for which we have 
no evolutionary tracks with rotation at subsolar 
chemical composition. In this case, our predicted 
RSG numbers are based on models with the old 
Geneva tracks only. This leaves open the possibil- 
ity that the RSG properties are indeed reproduced 
correctly in both SMC and LMC once the proper 
evolution models become available. 

Recently, iLevesque et al. 



new sample of RSGs in the SMC and LMC. Us- 
ing the same MARCS atmosphere models as in 



Levesque et al. ( 20051 ) they 



argue 



in favor of 

a 10% hotter effective temperature scale. Fur- 
thermore, metallicity-dependent systematic dif- 
ferences in the scale are found. Using this new ef- 
fective temperature scale, the observations match 
much better the theoretical predictions from the 
Geneva models as shown in Figs. l24l and l25l 

We have also compared the predictions for mod- 
els with solar chemical composition to observa- 
tions of Galactic RSGs with t he new effective 
tempe rature scale determined bv ILevesque et a~ 
(2005, Fig. 126]) . The definition of blue and red su- 



pergiants is also indicated in this figure. In this 
case, models with the new tracks reproduce the 
regime of Galactic RSGs very well. A simple con- 
clusion from this figure is that most RSGs are pro- 
duced by stars with initial masses lower than 30 
M Q , explaining the mismatch at lower Z in the 
SMC and LMC cases. Fig. [26] supports the ear- 
lier suggestion that a revised effective tempera- 
ture scale in the Magellanic Clouds could help to 
solve the mismatch with the new models. How- 
ever, the state of art of the recent models with 
rotation is still undergoing revision, and improved 
models could reach the red part of the observa- 
tions. 

The complexity of the evolutionary models and 
the incomplete information about different as- 
pects of massive stars do not yet allow us to pro- 
duce a more complete set of evolutionary models. 



Maeder fc Mevnet I (|200lh produced a set of tracks 
of v mt = 300 km s _1 with Z = 0.004. These tracks 
were specifically built with the purpose to reach 
the RSG zone. In contrast, the tracks used in the 
present work were produced with the goal of repro- 
ducing the overall properties of all massive stars 
in the upper HRD. In particular, modeling the be- 
havior of WR stars was a major objective. In order 
to simplify the calculations, convection in the en- 
velope was performed using a density scale height 
instead of a pressure scale height. This choice af- 
fects the temperatures in the RSG phase, making 
the stars hotter. The very good agreement with 
Galactic RSG favors this choice at solar chemical 
composition. 



(2006) analyzed 
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4.3. Observable star number ratios 

Reproducing the duration of individual stel- 
lar phases, such as the O, WR, or RSG phase, 
i s a major quest in stel l ar ev olution modeling 
(jMevnet fc Maeder I liuM 120051 ). In this section 
we use the predicted ratios of Section 3. We as- 
sume that the star numbers are from a stable pop- 
ulation, i.e., star formation has been ongoing for a 
sufficiently long interval to allow the relevant stel- 
lar phases to reach equilibrium. The input param- 
eters are the same as those described in Section 3. 

4.3.1. WC/WN 

Fig. [23 shows the WC/WN ratio for di ffer- 



ent stellar systems from iMassev et al. I (|2003l ) for 



the SMC. IBrevsacher et al. I (|l999h for the LMC, 



Hadfield et al. I (12007 ) for the sola r neighborhood 
Crockett et all (|2006t ) for M 3 1. lHadfield et al. 



(2005) for M 83, and for IC 10 bv lMassev fc Holmes 
([2002) . We have included in the figure three mod- 
els with different set of tracks: old Geneva (solid 
line), new tracks of v rot = 300 km s _1 (dotted 
and dot-dashed lines), and new tracks of v Iot = 
km s _1 (short-dashed line). For these models we 
have adopted a temperature limit to enter to the 
WR phase of logT e ff = 4.0, and we have used 
a temperature limit of logT e g = 4.4 just for the 
new tracks of i> ro t = 300 km s _1 (dot-short dashed 
line). 

As shown in Section 3, very high values for 
WC/WN are obtained for models using the old 
Geneva tracks and a higher temperature limit of 
log T c ff = 4.4 for the WR phase. The previous 
values for this ratio are higher than the upper 
limit of the plot in Fig. [27l In contrast, the new 
tracks with rotation do no longer overestimate the 
WC/WN ratio at low Z. At high Z, however, 
they seem to underestimate this ratio. Note that 
the dotted line fits better the trend shown by the 
data, especially for low values of the oxygen abun- 
dance. 

Models with the new tracks of v Iot = 300 km 
s _1 predict a higher number of WN stars for so- 
lar and higher oxygen abundance than expected 
from the observations. As discussed in the pre- 
vious sections, the number of WNs has increased 
because (i) the most massive stars enter the WR 
phase earlier and (ii) this phase lasts longer due to 
a combination of rotation and mass loss (in partic- 



ular at higher metallicities where these effects are 
more important). Models with tracks of v m t = 
km s _1 predict a ratio in better agreement with 
the observation at solar composition. An alterna- 
tive way to re solve this apparent paradox has been 
suggested by IMevnet fc Maeder I (2005), who ar- 
gued in favor of an underestimate of the mass loss 
during the WNL phase. 

The WC/WN ratio for IC 10 is still unclear 
since the statistics was based on unconfirmed WR 
candidates, but seems to be in the range of the 
model prediction. The Milky Way data were based 
upon stars within galacto-centric distances < 7 
kpc, 7-10 kpc, and > 10 kpc. 

4.3.2. WR/O 

An important observational constraint for stel- 
lar evolution models is provided by the varia- 
tion of the number ratio of WR to O-type stars 
with metallicity. Fig. [28] shows the WR/O ra- 
tio for different stellar p opulations taken from 
Maeder fc Mevnet I (|l994h at different radius from 



the Sun. The value for the SMC agre es with recent 
resul ts with ~ 1000 O-type stars fevans et al~ 
|2004| ) versus 12 WR (jMassev et al. fl2003l) giving 
a WR/O = 0.01. The data are compared to mod- 
els with the old Geneva tracks (solid line), new 
tracks of v Iot = 300 km s _1 (dotted line), and new 
tracks of w ro t = km s _1 (short-dashed line). 

Symbols annotated with an N for the SMC and 
LMC are based on recent determinations of the 
oxygen a bundances for these t wo galaxies. For 
the SMC, ILennon et a"l~l (|2003h fou nd a value of 



12 + \og{0/H) = 8.15 in B stars. ICosta et al. 
(2000) derived an oxygen abundance of 8.22 in a 
sample of planetary nebulae. In this work we have 
adopted the oxygen abundance obtained when us- 
ing a value of 8.75 for the solar composition and 
scaling this value to Z = 0.004 (i.e., dividing by a 
factor of 5). We did the same for the LMC oxy- 
gen abundance. The fit to the observations by 
our models becomes slightly better if we adopted 
a lower solar oxygen abundance, i.e., shifting the 
zero point of the models along the horizontal axis. 

Recent determinations of the solar oxygen 
abundances give val ues significantly lower (12 



log(Q/g) = 8.69: lAllende Prieto et al. I l2002t 



Asplun. Grevesse. fc SauvaH 20051 ) than the one of 



Anders fc Grevesse ( 19891 ). The dot-long dashed 
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line shows the same values for the WR/O ratio 
but considering lower values for the oxygen abun- 
dance for solar and higher chemical composition 
in Fig. Hi 

Obviously, the real situation is more complex 
than implied by our simple scaling, in particu- 
lar since the Fe abundance is the more important 
parameter driving evolution. Nevertheless, there 
seems to be a lack of WR stars in the models for 
solar and higher chemical composition. 

4.3.3. R/WR 

Fig. shows our models for the ratio of RSGs 
over WR stars (R/WR) using the same layout as 
in Fig. [27| Sym bols are the observed ratios from 



Massey 



20021 ). The figure includes only RSGs 
that lie in the same areas for which surveys are 
complete for WRs and follow the definition of a 
RSG as described in Section 3.2.3. 

Models using the old Geneva tracks (solid line) 
do not match the observations, except for the 
LMC oxygen abundance. The new tracks of v mt = 
300 km s _1 do not include tracks for masses < 30 
Mq for oxygen abundance less than solar. The 
mass range above 30 M is too small to produce a 
significant number of RSGs with the new tracks. 
Therefore we generate the same theoretical values 
as those obtained with the old Geneva tracks. For 
this reason, we have truncated the graphs below 
solar oxygen abundance. 

For solar and higher oxygen abundances, the 
new models and the data disagree. This is the re- 
sult of too few WR stars being predicted by the 
models. The adopted temperature limit for the 
WR phase does not significantly affect the results. 
This can be seen as seen for the models with tracks 
of Urot = 300 km s^ 1 (dot-short dashed line). Re- 
cent results for M 33 and M 31 suggest lower oxy- 
gen abundances in HII regions by about .1 dex 



(|Massev et ai~l 120031: ICrockett et al. Il2006h . The 



lower abundances do not significantly affect our 
conclusions but should be kept in mind when in- 
specting Fig. [29] 

4.3.4. B/R 

The final star number ratio to compare to our 
models is the blue-to-red supergiant ratio. Fig. [30] 
shows the B/R ratio for different sets of data as 
a function of the oxygen abundance. There arc 



several determinations for this ratio in the Magel- 
lanic Clouds. We used the ho mogeneou s surve y of 
stars in the SMC and LMC by Massey I (|200l tri- 
angles) and iMassev fe Olsen I (|2003l . squares) for 
our comparison. 



Massev fe Olsen 



(|2003f) discuss contamination 
effects for K and M-type stars in surveys when a 
color cut-off of (V — R) > 0.6 is imposed. Fore- 
ground contamination can significantly affect the 
observed numbers. 11% and 5.3% of the bona 
fide RSGs seen towards the SMC and LMC, re- 
spectively, are foregroun d objects in the survey of 



Massey fe Olsen I (|2003l ). Another error source is 
the calibration of the relation to transform colors 
into effective temperature. 

We al so included in Fig [30l th e results of the 
study of lLanger fe Maeder I (|1995f ) who analyzed 
different samples of blue and red supergiants in as- 
sociations and in extragalactic fields. The trend of 
these observations suggests a lower value of B/R 
for lower chemical composition and higher values 
for higher chemical composit ion. Thi s resul t dis- 
agrees with the numbe rs of IMassev I (|2002l ) and 



Massev fe Olsen I (|2003f) 



The third data set i n Fig.JSDlhas been obtained 



by lEggenberger et aT~ ( 20021 ). Their fit was de- 
rived by using a sample of 45 clusters in the Milky 
Way and 4 clusters in Magellanic Clouds. To 
obtain the long-dashed line in the plot, we as- 
sumed a value of B/R = 30 at solar composi- 
tion, which is simi l ar to the value obtained by 
Eggenberger et al~l ( 20021 ) for Galactic clusters. 

As before, three models are shown in Fig. [30] 
with the line types as in Fig. [28] We do not plot 
the B /R ratio for chemical composition lower than 
solar for the same reason as in the previous sub- 
section. Models with the old Geneva tracks do 
not fit the observations, whereas models with the 
new tracks of v ro t = 300 km s _1 match the obser- 
vations for oxygen abundance higher than solar. 
Models with the new tracks of u ro t = km s _1 
predict similar values as the models with the old 
Geneva tracks. As argued before, this ratio does 
not depend on the minimum temperature for the 
WR phase. 

It is evident that the data are very inconsis- 
tent for subsolar chemical composition. Neverthe- 
less, the data agree in a trend of higher B/R at 
higher O/H. Our models predict a similar trend. 



15 



However, our predictions are still affected by un- 
certainties in the zero point for the oxygen abun- 
dance. As in the case of the WR/O ratio, our val- 
ues agree better with the observations if we adopt 
a lower solar oxygen abundance as shown by the 
dot-long dashed line in Fig. [30] 

5. Conclusions 

Progress in stellar evolution modeling is a 
rather discontinuous process when viewed from 
a "user's" perspective. Population synthesis mod- 
elers benefit from and are at the mercy of updates 
and new releases by the major stellar evolution 
modeling groups. Evolution models for massive 
stars have undergone dramatic improvements over 
the past several decades. Major milestones were 
set with the realization of the importance of stel- 
lar mass loss, of convective processes, and of the 
influence of chemical composition. Each time the 
release of a new set of stellar evolution models led 
to sometimes major revisions of the predictions by 
population synthesis models. 

The latest major release of revised evolution 
models by the Geneva and Padova groups occurred 
in the 1990s when new opacities were included and 
the importance metallicity-dependent mass loss 
was highlighted. These sets of evolution models 
are currently employed by virtually all population 
synthesis code in use. Now we are in an epoch 
of another quantum leap in massive star evolu- 
tion. The importance of stellar rotation for mas- 
sive stars has been recognized and evolution calcu- 
lations have been done over the past years. Only 
very recently has a grid of evolution models with 
rotation been completed that is sufficiently broad 
and dense to be suitable for implementation. This 
model grid, with ZAMS velocities of and 300 km 
s _1 , has been implemented into Starburst99. 

We acknowledge the shortcomings of the mod- 
els presented in this paper. The limited mass cov- 
erage of the available stellar models often intro- 
duces serious biases making the comparison with 
observations difficult. Nevertheless, some of the 
trends we have found in our study are significant 
improvements to the stellar populations study. 

Massive hot, blue stars tend to be more lumi- 
nous and hotter when rotation is accounted for. 
This results from the increased convective core 
and larger amount of mixing reducing the sur- 



face opacities. When coupled with synthesis mod- 
els, these properties lead to significantly different 
population parameters. For a fixed stellar mass, 
a young population with rotating stars becomes 
more luminous at almost all wavelengths of inter- 
est. Specifically, the ultraviolet to near-infrared 
luminosities can increase by factors of order 2, in 
particular when RSGs are present. If such models 
are compared with observed light-to-mass ratios of 
young star clusters, the derived mass may become 
smaller by the same factor in comparison with the 
results of earlier studies. This has important con- 
sequences for our understanding of the low-mass 
end of the stellar initial mass function. 

The impact of the new models at wavelengths 
below 912 A is even more dramatic. The SED 
in the extreme ultraviolet becomes enhanced and 
harder when compared with previous synthesis 
models. The effect is most pronounced in the ion- 
ized helium continuum below 228 A due to the 
increased WR star contribution. Even at shorter 
wavelengths in the neutral helium and hydrogen 
continua below 504 and 912 A, respectively, the 
changes are significant. Since the Lyman contin- 
uum flux increases by a factor of 2 to 3, the de- 
rived stellar masses or star-formation rates from 
hydrogen recombination will decrease by the same 
factor. A large fraction of the increased photon 
output is provided by WR stars whose frequency 
is strongly dependent on the chemical composi- 
tion. Therefore the prediction recombination-line 
fluxes and star-formation tracers become more de- 
pendent on Z than in previous models. 

Given the potential impact of these new predic- 
tions it becomes crucial to perform a critical as- 
sessment of the additional assumptions that were 
made. The increase and hardening of the spec- 
trum is the result of the stellar interior and atmo- 
sphere connection, which rests on the definition of 
the effective temperature and radius of the star. 
Moreover, atmosphere models for WR stars must 
be adopted. Another major area of concern are 
the mass-loss rates, in particular during the WR 
phase. The mass loss rates used in the default ver- 
sion of the new tracks with rotation do not include 
a metallicity dependence of the rates. It is also a 
matter of debate how the competing effects and 
mass loss affect the stellar properties. 

Nevertheless, this study should be viewed as a 
first step in the development of more realistic and 
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accurate models for stellar populations. Should 
our initial results be confirmed when more mature 
evolutionary tracks are available, several widely 
used star-formation and initial mass function in- 
dicators will have to be revised. 

Thanks to Leonardo Ubeda who helped with 
the Hess diagrams. Support for this work was pro- 
vided by NASA through grant 01166 from the 
Space Telescope Science Institute, which is op- 
erated by the Association of Universities for Re- 
search in Astronomy, Inc., under NASA contract 
NAS5-26555. 
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Table 1 

Minimum mass limits for entering the WR phase for different sets of tracks (in 



Chemical Composition 



Mass limit 


Z = 0.004 


Z = 0.008 


Z - 0.02 


Z = 0.04 


M WR a (u rot = 300 km s _1 ) 


32 


25 


22 


21 


MwR(»rot = km s _1 ) 


b 




37 


29 


Mowh" (v [0 t = 300 kins" 1 ) 


75 


69 


45 


39 


MowR(»rot = km s _1 ) 






62 


42 


M WR (HML C ) 


42 


35 


25 


21 


M W r(LML c ) 


52 


42 


32 


25 



a MowR- The mass limit to enter to the WR phase on the main-sequence. Mwr. 
The mass limit to enter the WR phase at any time during the evolution of the star 
IMevnet fc Maeder I [20031 [20051) . 

b No computations arc available with v TO t — km s — 1 

c 01d Geneva tracks denoted as HML and LML for high and normal mass-loss rates, 
respectively. 



Table 2 

Mass range (in M q ) for the new models with rotation. 



Tracks with initial v ro t — 300 km s 



z 


= 0.004 Z = 0.008 Z = 0.02 


z = 


0.04 




120 120 120 




120 




85 




85 




60 60 60 




60 




40 40 40 




40 




30 30 








25 




25 




20 




20 




15 








12 








9 








Tracks with initial v vo t — km 


s- 1 




z 


= 0.004 Z = 0.008 Z = 0.02 


z = 


0.04 




120 




120 




85 








60 




60 




40 








25 




25 




20 








15 








12 








9 








Tracks with mctallicity dependent WR rates 


z 


= 0.004 Z = 0.008 Z = 0.02 


z = 


0.04 



120 
85 
60 
40 



21 



Table 3 

smc observations for comparison with z — 0.004 models. 



Name 


logT cff 


log L/L B 


12 + logN/H 


V sin i 


Mass 


dM/dt 












[km s" 1 ] 


[M Q ] 


[10- 6 M o yr 








With UVES at VLT Trundle et a 


. (2004) 








AV215 


4.431 


5.63 ± 0.06 


7.94 ± 0.10 




90 


26 


1.3500 ± 


14 


AVI 04 


4.439 


5.31 ± 0.06 


7.40 ± 0.13 




70 


19 


0.4000 ± 


05 


AV216 


4.415 


5.00 ± 0.06 


7.63 ± 0.19 




75 


36 


0.1200 ± 


02 


AV362 


4.146 


5.50 ± 0.06 


8.20 ± 0.25 




50 


17 


0.8000 ± 


08 


AV22 


4.161 


5.04 ± 0.06 


7.90 ± 0.14 




40 


8 


0.2300 ± 


02 


AV18 


4.279 


5.44 ± 0.06 


7.50 ± 0.31 




40 


17 


0.2300 ± 


03 


AV210 


4.312 


5.41 ± 0.06 


7.60 ± 0.12 




65 


15 


0.2000 ± 


05 


Skl91 


4.352 


5.77 ± 0.06 


7.63 ± 0.20 




95 


33 


0.6800 ± 


07 






EMMI at the 


NTT Trundle & Lennon (2005 


) 






AV78 


4.332 


5.92 ± 0.06 


8.30 ± 0.26 






57 


2.2900 ± 


34 


AV264 


4.352 


5.44 ± 0.06 


7.88 ± 0.12 






16 


2900 + 


06 


AV242 


4.398 


5.67 ± 0.06 


7.39 ± 0.19 






35 


0.8400 ± 


13 


AV420 


4.431 


5.35 ± 0.06 


7.44 ± 0.12 






19 


0.3400 ± 


15 


AV151 


4.204 


5.28 ± 0.06 


7.55 ± 0.18 






15 


0.1600 ± 


04 


AV56 


4.217 


5.88 ± 0.06 


7.93 ± 0.22 






38 


0.5100 ± 


08 


AV443 


4.217 


5.79 ± 0.06 


8.26 ± 0.33 






30 


0.4500 ± 


09 


AV10 


4.230 


5.21 ± 0.06 


7.69 ± 0.27 






13 


0.1500 ± 


03 


AV373 


4.279 


5.42 ± 0.06 


7.46 ± 0.00 






16 


0.1600 ± 


04 


AV96 


4.342 


5.39 ± 0.06 


7.69 ± 0.24 






15 


0.2400 ± 


06 






O-tvpc dwarfs from Bouret et al. 


(2003) 








MPG 355 


4.72 


6.03 


7.95 




110 


65 


2.5100 




MPG 324 


4.62 


5.52 


7.27 




70 


40 


0.2500 




MPG 368 


4.60 


5.42 


7.75 




60 


40 


0.1600 




MPG 113 


4.60 


5.13 


7.27 




36 


30 


0.0030 




MPG 487 


4.54 


5.13 


6.41 




20 


25 


0.0030 




MPG 12 


4.48 


4.93 


7.95 




60? 


20 


0.0001 








O-typc supcrgiants from Hillicr et al. (2003 










4.53 


5.60 


6.30 




60 










4.52 


5.53 


8.39 




70 










AF-type 


supcrgiants from Venn (1999) 


and Venn & Przvbilla 


(2003) 






3.98 


4.77 


7.79 




25 










3.975 


4.59 


7.70 




15 










3.93 


4.47 






25 










3.92 


5.18 


7.72 




20 










3.91 


4.85 


7.88 




15 










3.905 


4.62 


7.07 




25 










3.905 


4.87 


7.29 




20 










3.895 


4.98 


7.56 




20 










3.875 


4.58 


7.80 




20 










3.86 


4.48 


6.80 




36 
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Table 4 



Average deprojected rotational veloci ties d erived from the data of IConti fc Ebbets 

(|l977h . 



L/L 



4.445 


4 


460 


97 


4.485 


5 


630 


104 


4.525 


5 


435 


107 


4.505 


5 


340 


108 


4.420 


4 


815 


108 


4.415 


5 


415 


112 


4.535 


5 


790 


112 


4.515 


5 


725 


116 


4.465 


5 


570 


119 


4.570 


5 


905 


131 


4.610 


5 


905 


134 


4.555 


5 


585 


144 


4.580 


5 


730 


153 


4.645 


6 


205 


153 


4.525 


5 


760 


159 


4.630 


5 


760 


164 


4.680 


5 


995 


166 


4.465 


5 


040 


166 


4.540 


5 


530 


167 


4.500 


4 


850 


175 


4.535 


5 


110 


175 


4.485 


5 


190 


179 


4.650 


5 


900 


180 


4.485 


4 


715 


181 


4.545 


5 


820 


193 


4.615 


5 


625 


208 


4.580 


5 


415 


221 


4.570 


5 


645 


22!) 


4.600 


5 


815 


238 


4.555 


5 


230 


243 


4.520 


4 


985 


251 


4.600 


5 


520 


268 


4.570 


5 


325 


270 
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Table 5 

Same as in Table 0] but derived from the data of IHuang fc Gies I (120061 ) . 



T ft 


L/L e 


V ro t 


Teff 


L/L e 


^rot 




L/L e 


V ra t 


4.400 


4 


500 


1 


4 


.250 


2 . 


750 


131 


4 


.175 


2 


.500 


192 


4.350 


4 


250 


6 


4 


.150 


3 


000 


132 


4 


.200 


2 


.750 


193 


4.375 


4 


.750 


8 


4 


.250 


3. 


000 


136 


4. 


.400 


3 


.875 


199 


4.275 


4. 


000 


14 


4 


.350 


3. 


625 


136 


4 


.325 


3 


.625 


202 


4.325 


4 


.375 


23 


4 


.450 


4 


.375 


130 


4 


.250 


3 


.250 


203 


4.200 


3 


250 


20 


4 


.200 


2. 


625 


140 


4 


.450 


4 


.500 


205 


4.300 


3. 


000 


31 


4 


.475 


4. 


125 


144 


4 


.175 


2 


.750 


208 


4.225 


2 


.750 


36 


4 


.100 


2 


.875 


144 


4 


.375 


3 


.875 


213 


4.150 


3 


250 


39 


4 


.325 


3. 


500 


145 


4 


.250 


3 


.125 


213 


4.450 


4 


625 


41 


4 


.450 


4. 


125 


146 


4 


.375 


4 


.125 


214 


4.225 


3 


.375 


45 


4 


.350 


3. 


.125 


140 


4 


.250 


3 


.375 


218 


4.050 


2 


.750 


47 


4 


.325 


3 


250 


140 


4 


.225 


3 


.125 


222 


4.375 


3. 


.500 


47 


4 


.400 


4. 


250 


151 


4 


.275 


3 


.250 


223 


4.425 


4. 


.375 


52 


4 


.150 


2. 


500 


153 


4 


.225 


3 


.500 


228 


4.300 


4 


000 


60 


4 


.425 


3 


.875 


155 


4 


.225 


2 


.875 


230 


4.150 


2 


.750 


62 


4 


.350 


4. 


375 


156 


4 


.475 


4 


.750 


232 


4.125 


2 


500 


66 


4 


.125 


2 


.875 


150 


4 


.275 


3 


.375 


236 


4.250 


2 


625 


74 


4 


.200 


3 


000 


150 


4 


.225 


3 


.000 


240 


4.425 


3 


.750 


78 


4 


.400 


4 


.125 


150 


4 


.325 


3 


.750 


253 


4.300 


3. 


250 


86 


4 


.275 


3 


000 


162 


4 


.325 


4 


.000 


255 


4.250 


3 


.500 


88 


4 


.350 


3. 


875 


162 


4 


.425 


4 


.625 


255 


4.375 


3 


625 


91 


4 


.275 


3 


.125 


163 


4 


.275 


3 


.500 


263 


4.200 


3. 


.125 


93 


4 


.425 


4. 


250 


167 


4 


.175 


2 


.625 


269 


4.425 


4 


500 


93 


4 


.350 


3. 


750 


167 


4 


.225 


2 


.625 


271 


4.425 


4 


000 


94 


4 


.400 


3 


750 


170 


4 


.200 


2 


.500 


271 


4.175 


3 


000 


95 


4 


.050 


2. 


625 


171 


4. 


.325 


3 


.875 


277 


4.300 


3 


.375 


102 


4 


.300 


3. 


625 


171 


4 


.350 


4 


.000 


278 


4.350 


3. 


.375 


103 


4 


.425 


4 


.125 


171 


4 


.300 


3 


.750 


281 


4.325 


4 


250 


106 


4 


.375 


3 


750 


172 


4 


.275 


3 


.625 


292 


4.475 


4 


.875 


110 


4 


.350 


3 


500 


172 


4 


.375 


4 


.000 


293 


4.150 


3 


500 


118 


4 


.400 


4 


.375 


172 


4 


.325 


4 


.500 


297 


4.450 


4 


000 


122 


4 


.300 


3. 


500 


172 


4 


.200 


2 


.875 


309 


4.400 


4 


000 


126 


4 


.250 


2 


875 


173 


4 


.275 


3 


.750 


326 


4.350 


4 


.125 


127 


4 


.475 


4 


625 


180 


4. 


.050 


2 


.500 


351 


4.450 


4 


250 


120 


4 


.375 


4. 


.375 


183 


4 


.075 


3 


.000 


397 


4.150 


2. 


.625 


130 


4 


.325 


3 


.375 


180 


4 


.075 


2 


.500 


402 
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4.8 4.6 4.4 4.2 4 3.8 4.8 4.6 4.4 4.2 4 3.8 

'°e(T.rr> l°e(T,„) 



Fig. 2. — Isochrones produced with Starburst99. 
Z = 0.008: ages of 0.01, 1, 2, 3, 4, 5, 6, and 7 Myr. 
Z = 0.020: ages of 0.01, 1, 2, 5, 8, 11, 14, 17, 
20, 23, 26, and 29 Myr. The oldest isochrone in 
each figure corresponds to the lifetime of the least 
massive star included in the tracks with v TO t — 
300 km s" 1 . 



Fig. 1. — Comparison of three 60 Mq tracks with 
different physics input from the Geneva group; 
solid: the original high mass loss, short-dashed: 
Vrot = 0, dotted: u rot = 300 km s _1 . The triangle 
indicates the beginning of the WR phase for the 
u ro t = 300 km s _1 models. 
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Fig. 3. — SEDs for young populations with solar 
composition and 10 6 M Q . Solid: previous tracks 
with enhanced mass loss; short-dashed: v rot = 0; 
dotted: v mt = 300 km s~ . For clarity, the SEDs 
for 3, 4, 5, and 6 Myr are displaced by -10, -20, -30 
and -40 dex from the original value, respectively. 
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Fig. 5.— Red part of the SEDs for ages 11 - 19 
Myr (top to bottom) and displaced incrementally 
by -1.0 dex from the original luminosity, except for 
the 11 Myr old model. Line types as in Fig.[3l 



Fig. 4. — Same as Fig. [3J but for ages of 7, 8, 9, 
10, and 15 Myr, displaced by -10, -20, -30, and -40 
dex, respectively, except for the model at 7 Myr 
which is mishitted. 
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Fig. 6. — Ratio of WR over O stars for an instan- 
taneous burst. Same line types as in the previ- 
ous figures: logT G ff > 4.4 (thick solid and dotted 
lines) and logT c ff > 4.0 (long-dashed, dot-dashed 
and thick short-dashed lines) for models with the 
old Geneva tracks, new tracks with v Iot = 300 
and v m t = km s , respectively. Thick black 
dot-long-dashed line: models with new tracks for 
different WR mass-loss rates and log T c g > 4.0 at 
Z = 0.040. 
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Fig. 7. — Same as Fig. [SJ but for continuous star 
formation. 
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Fig. 8. — Same as Fig. [6l but for the ratio of WC 
over WN stars. 
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Fig. 9. — Same as Fig. [7l but for the ratio of WC 
over WN stars. 
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r ] r 1 log(Age [Myr]) log(Age [Myr]) 

0.5 1 1.5 2 0.5 1 1.5 2 

log(Age [Myr]) log(Age [Myr]) 

Fig. 13.— Afboj, Liboo, M B , M v , Mj, and M K vs. 

Fig. 11.- Same as Fig. but for the ratio of time for models with the old Geneva ( red ' solid )' 

blue over red supergiants. Note that the apparent new Geneva wlth v ^ = ( black ' short-dashed), 

increase of B/R with decreasing Z seen for Z = and new models wlth v ^ = 300 km s ( blue ' dot " 

0.004 and 0.008 is an artifact of incomplete tracks. ted ) tracks ' The bands most affected b y the new 

models are those in the ultraviolet and infrared. 

Solar chemical composition. 
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Fig. 14. — Color evolution for selected niters. 
Lines are as in Fig. Q21 a) Z = 0.004; b) Z = 
0.008; c) Z = 0.020. Model with Z = 0.020 and 
the new tracks of u ro t = 300 km s _1 evolve un- 
til 100 Myr whereas the other models terminate 
earlier. This is due to the larger mass range cov- 
ered (9 - 120 M Q ) than models at lower chemical 
composition. 
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Fig. 15. — Evolution of the mass-loss rate for 
five sets of models with 10 6 M . Four different 
chemical compositions are plotted. The figure for 
Z = 0.040 includes models with a mass-loss rate 
oc Z - 5 for the WR phase. Close to the ZAMS, 
the mass-loss rate is briefly higher in models with 
Vrot = 300 km s _1 (blue, dotted) than in models 
with the old tracks (red, solid). The short-dashed 
lines are models with the new tracks of u rot = km 
s _1 . The adopted stellar mass- loss rates for the ro- 
tating and non-rotating models are the same. A 
correction factor has been applied to the rotating 
models as explained in Sect. 2.1. 



logiAge [Myr]) log(Age [Myr]) 

Fig. 16. — Evolution of the mass loss for the same 
models as in Fig. [15] The total mass loss is similar 
in all models. 




0.5 i 1,5 
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Fig. 17. — The light-to-mass ratio for Lboi, £v> 
and Lk for different chemical compositions, a) 
Z = 0.020, and b) Z = 0.040. Same line types as 
before. Models with the new tracks of i> ro t = 300 
km s _1 are typically more luminous than those 
with the old Geneva tracks. 
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Fig. 19. — Same as in Fig |T8l but for the current 
mass. 




log[(N/H)/(N/H) ini }. The crosses indicate i°g(T.,,M) 

the validity limit of the models with the new 

tracks of v Iot = 300 km s -1 . Z = 0.004. Symbols Fig. 20. — Same as in FigHH but for the mass-loss 

denote data. rate. 
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Fig. 25.— Same as in Fig. El but for LMC RSGs 
and models with Z = 0.008. 
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Fig. 27.— Ratio of WC/WN as a function of 
the oxygen abundance. Models with a tempera- 
ture limit of logTcff = 4.0 for entering the WR 
phase and the old Geneva tracks (solid), new 
tracks of v Iot = 300 km s _1 (dotted), and new 
tracks of v ro t — km s _1 (short-dashed) are 
shown. An additional model with a tempera- 
ture limit of log T c ff = 4.4 and with the new 
tracks of v ro t = 300 km s _1 (dot-short dashed) 
is also shown. The data for the SMC, IC 10, 
LMC, M 33 and M 3 1 MW, and M 83 a re from 
Massev et al. I (120031): iMassev fc Holmes 



Brcysachcr ct al. (19991); ICrockett et al 



(2002); 
( 2006 ): 

Hadfield et all (|2007l) and lHadfield et al. I (|2005l ). 



00 3.75 3.50 



respectively. 



I°9(T,„ [K]) 



Fig. 26. — Same as in Fig. [25] but for Galactic 
RSGs and models with Z = 0.020. In this case, 



the T e ff and luminosity l imits are from IMassev 
(2002) and the data from lLevesaue et al. ( 20051 ). 
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Fig. 28. — Ratio of WR/O as a function of the 
oxygen abundance. Models are as in Fig. [27] but 
we have omitted the model with a higher tem- 
perature limit for entering the WR phase. The 
dot-long dashed line shows the trend of the mod- 
els with the new tracks of v Iot = 300 km s" 1 if 
the oxygen abundance in the models is decreased. 
The lower oxygen abundance is closer to the ob- 
servations. 




12+log(0/H) 



Fig. 29.— Ratio of RSGs over WR stars Model s 



are as in Fig. [27] Data are from iMassev I ([2003]). 
The new tracks with v Iot = 300 km s _1 for chem- 
ical composition lower than solar do not reach 
RSGs as defined by Massey. For solar and super- 
solar chemical composition, the models produce 
too few WR stars to match the observations. 
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Fig. 
giants. 

are data from | Massev 



30. — Ratio 
Models 



Massev fc Olsenl tl2003h 



tions of lLanger fc M aedcr 
stars of lLanger fc Maeder 



of blue over red super- 
are as in Fig. [28] Triangles 
squares are from 
ons are associa- 
circles are field 
The 



(12002), 
>ent as 



Eientago 
(1995). 



(1995). 



long- 
dashed line is a fit bv lEggenberger et al. I (|2002l ) to 



a large sample of clusters in different system s in- 
cluding those from[Langer^^aeder] (1995[). The 
dot-long dashed line denotes models with the new 
tracks of v rot — 300 km s _1 using updated values 
of the solar oxygen abundance. 
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